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Abstract. We review the present status of three-flavour neutrino oscillations, taking 
into account the latest available neutrino oscillation data presented at the Neutrino 
2008 Conference. This includes the data released this summer by the MINOS 
collaboration, the data of the neutral current counter phase of the SNO solar neutrino 
experiment, as well as the latest KamLAND and Borexino data. We give the updated 
determinations of the leading 'solar' and 'atmospheric' oscillation parameters. We find 
from global data that the mixing angle 9i3 is consistent with zero within 0.9(t and we 
derive an upper bound of sin^ 6*13 < 0.035 (0.056) at 90% CL {3a). 
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1. Introduction 



Thanks to the synergy amongst a variety of experiments involving solar and atmospheric 
neutrinos, as well as man-made neutrinos at nuclear power plants and accelerators jl| 
neutrino physics has undergone a revolution over the last decade or so. The long-sought- 
for phenomenon of neutrino oscillations has been finally established, demonstrating that 
neutrino flavor states (z/g, z/^, z/^) are indeed quantum superpositions of states (z^i, 1^2, 1^3) 
with definite masses The simplest unitary form of the lepton mixing matrix 

relating flavor and mass eigenstate neutrinos is given in terms of three mixing angles 
(^12, 6*13, 6*23) and three CP-violating phases, only one of which affects (conventional) 
neutrino oscillations j^. Here we consider only the effect of the mixing angles in 
current oscillation experiments, the sensitivity to CP violation effects remains an open 
challenge for future experiments 0, Together with the mass splitting parameters 
Am2i = ^2 — ml and Amg^ = r/ig — ml the angles 612, 623 are rather well determined by 
the oscillation data. In contrast, so far only upper bounds can be placed upon 613, mainly 
following from the null results of the short-baseline CHOOZ reactor experiment [gI with 
some effect also from solar and KamLAND data, especially at low ^m"^^ values [7|]. 

Here we present an update of the three-flavour oscillation analyses of Refs. 
and This new analysis includes the data released this summer by the MINOS 
collaboration J9|, the data from the neutral current counter phase of the SNO experiment 
(SNO-NCD) [lOj, the latest KamLAND ll| and Borexino 12|] data, as well as the results 
of a recent re- analysis of the Gallex/GNO solar neutrino data presented at the Neutrino 
2008 conference 
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In Section [2] we discuss the status of the parameters relevant for 
the leading oscillation modes in solar and atmospheric neutrinos. In Section |3] we present 
the updated limits on 6*13 and discuss the recent claims for possible hints in favour of a 
non-zero value made in Refs. [13,0, 16]. We summarize in Section |H 



2. The leading 'solar' and 'atmospheric' oscillation parameters 

Let us first discuss the status of the solar pararneters 612 and Am|p The latest data 



release from the KamLAND reactor experiment Uj has increased the exposure almost 
fourfold over previous results [13] to 2.44 x 10'^^ proton-yr due to longer lifetime and an 
enlarged fiducial volume, corresponding to a total exposure of 2881 ton-yr. Apart from 
the increased statistics also systematic uncertainties have been improved: Thanks to 
the full volume calibration the error on the fiducial mass has been reduced from 4.7% 
to 1.8%. Details of our KamLAND analysis are described in appendix A of Ref. 0]. We 
use the data binned in equal bins in l/E" to make optimal use of spectral information, 
we take into account the (small) matter effect and carefully include various systematics 
following Ref. 
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As previously, we restrict the analysis to the prompt energy range 
above 2.6 MeV where the contributions from geo-neutrinos 191] as well as backgrounds 
are small and the selection efficiency is roughly constant HI]. In that energy range 1549 



reactor neutrinos events and a background of 63 events are expected without oscillations. 
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whereas the observed number of events is 985 20 1. 

The Sudbury Neutrino Observatory (SNO) has released the data of its last phase, 
where the neutrons produced in the neutrino NC interaction with deuterium are detected 
mainly by an array of '^He proportional counters to measure the rate of neutral- 
current interactions in heavy water and precisely determine the total active boron 



solar neutrino flux, yielding the result 5.54^Q;3^(stat)i'iQ'3^(syst) x 10^cm~^s~^ 10|. The 
independent ^He neutral current detectors (NCD) provide a measurement of the neutral 
current flux uncorrelated with the charged current rate from solar z/g, different from the 
statistical CC/NC separation of previous SNO phases. Since the total NC rate receives 
contributions from the NCD as well as from the PMTs (as previously) a small (anti-) 
correlation between CC and NC remains. Following Ref. 16| we assume a correlation of 
p = —0.15. In our SNO analysis we add the new data on the CC and NC fluxes to the 
previous results 21| assuming no correlation between the NCD phase and the previous 
salt phase, see Ref. 0] for further details. The main impact of the new SNO data is due 
to the lower value for the observed CC/NC ratio, (0cc/0nc}1^° = 0.301 ± 0.033 jlO 



compared to the previous value {4>cc/(pNcT^^^ = 0.34 ±0.038 [2l|. Since for neutrinos 
0cc/0NC ~ Pee ~ siu^ 6*12, adding the new data point on this ratio with the lower value 
leads to a stronger upper bound on sin^6'i2. 

We also include the direct measurement of the ^Be solar neutrino signal rate 
performed by the Borexino collaboration 12|. They report an interaction rate of 
the 0.862 MeV ^Be neutrinos of 49±3(stat)±4(syst) counts/(day-100 ton). This 
measurement constitutes the first direct determination of the survival probability 
for solar z/g in the transition region between matter-enhanced and vacuum-driven 
oscillations. The survival probability of 0.862 MeV ^Be neutrinos is determined to 
be P^^^'"^"^ = 0.56 ± 0.1. We find that with present errors Borexino plays no significant 
role in the determination of neutrino oscillation parameters. Apart from the fact that 
the uncertainty on the survival probability is about a factor 3 larger than e.g., the 
uncertainty on the SNO CC/NC ratio measurement, it turns out that the observed 
value for Pee quoted above practically coincides (within 0.1 a) with the prediction at the 
best fit point: pjBe,pred ^ q_55_ 

The new data from SNO and Borexino are combined with the global data on solar 



neutrinos 



22|, |23|, 124, l25| , where we take into account the results of a recent re- analysis 
of the Gallex data yielding a combined Gallex and GNO rate of 67.6±4.0±3.2 SNU 13 . 

Fig. [1] illustrates how the determination of the leading solar oscillation parameters 
6i2 and Amgi emerges from the complementarity of solar and reactor neutrinos. From 
the global three-flavour analysis we find (la errors) 

sin^ Ou = 0.304+°:°22 ^ ^^2^ _ 7.65l[^:| x 10"^ eV^ . (1) 

The numerical profiles shown in Fig. [1] have to very good accuracy the Gaussian 
shape = ~ a^best)^/c"^, when the different a for upper an lower branches are 
used as given in Eq. ([1]). Spectral information from KamLAND data leads to an 
accurate determination of Amgi with the remarkable precision of 8% at 3a, defined 
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Figure 1. Determination of the leading "solar" oscillation parameters from the 
interplay of data from artificial and natural neutrino sources. We show x^-profiles 
and allowed regions at 90% and 99.73% CL (2 dof) for solar and KamLAND, as well 
as the 99.73% CL region for the combined analysis. The dot, star and diamond indicate 
the best fit points of solar data, KamLAND and global data, respectively. We minimise 
with respect to Am^i, 623 and ^13, including always atmospheric, MINOS, K2K and 
CHOOZ data. 



as (x^PP^"' - a;'°™^'')/(x"PP^'' + We find that the main hmitation for the Am^^ 

measurement comes from the uncertainty on the energy scale in KamLAND of 1.5%. 
KamLAND data start also to contribute to the lower bound on sin^6'i2, whereas the 
upper bound is dominated by solar data, most importantly by the CC/NC solar neutrino 
rate measured by SNO. The SNO-NCD measurement reduces the 3a upper bound on 
sin^^ia from 0.40 B to 0.37. 

Let us now move to the discussion of the status of the leading atmospheric 
parameters 623 and Amg^^. The Main Injector Neutrino Oscillation Search experiment 
(MINOS) has reported new results on z/^ disappearance with a baseline of 735 km based 
on a two-year exposure from the Fermilab NuMI beam. Their data, recorded between 
May 2005 and July 2007 correspond to a total of 3.36 x lO^o protons on target (POT) Q, 



more than doubling the POT with respect to MINOS run I 26|, and increasing the 




Figure 2. Determination of tlie leading "atmospheric" oscillation parameters from 
the interplay of data from artificial and natural neutrino sources. We show x^-profiles 
and allowed regions at 90% and 99.73% CL (2 dof) for atmospheric and MINOS, as 
well as the 99.73% CL region for the combined analysis (including also K2K). The dot, 
star and diamond indicate the best fit points of atmospheric data, MINOS and global 
data, respectively. We minimise with respect to Am'^i, O12 and ^13, including always 
solar, KamLAND, and CHOOZ data. 



exposure used in the latest version of Ref. |8[ by about 34%. The latest data confirm the 
energy dependent disappearance of z/^, showing significantly less events than expected 
in the case of no oscillations in the energy range < 6 GeV, whereas the high energy part 
of the spectrum is consistent with the no oscillation expectation. We include this result 
in our analysis by fitting the event spectrum given in Fig. 2 of Ref. jol. Current MINOS 
data largely supersedes the pioneering K2K measurement 27| which by now gives only 
a very minor contribution to the Am|^ measurement. 

We combine the long-baseline accelerator data with atmospheric neutrino 
measurements from Super-Kamiokande 28|, using the results of Ref. [Sj, see references 
therein for details. In this analysis sub-leading effects of Am^i in atmospheric data are 



neglected, but effects of 6*13 are included, in a similar spirit as in Ref. [29 



Fig. [2] illustrates how the determination of the leading atmospheric oscillation 
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Figure 3. Constraints on sin 613 from different parts of the global data. 



parameters 62^ and |Am|]^| emerges from the complementarity of atmospheric and 
accelerator neutrino data. We find the following best fit points and la errors: 



sin^ e 



23 



0.50: 



-0.07 
-0.06 ' 



\Aml^\ 



'eV' 



(2) 



The determination of lAmg^^l is dominated by spectral data from the MINOS long- 
baseline z/^ disappearance experiment, where the sign of Amg^ (i.e., the neutrino mass 
hierarchy) is undetermined by present data. The measurement of the mixing angle 623 
is still largely dominated by atmospheric neutrino data from Super-Kamiokande with a 
best fit point at maximal mixing. Small deviations from maximal mixing due to sub- 



leading three-flavour effects have been found in Refs. [30|, l3l| , see, however, also Ref . [32 



for a preliminary analysis of Super-Kamiokande. A comparison of these subtle effects 
can be found in Ref. 33|]. At present deviations from maximality are not statistically 
significant. 



3. Status of 

The third mixing angle ^13 would characterize the magnitude of CP violation in neutrino 
oscillations. Together with the determination of the neutrino mass spectrum hierarchy 
(i.e., the sign of Am'^i) it constitutes a major open challenge for any future investigation 
of neutrino oscillations 0, ^ . 

Fig. [3] summarizes the information on 613 from present data. Similar to the case of 
the leading oscillation parameters, also the bound on 6*13 emerges from an interplay of 
different data sets, as illustrated in the left panel of Fig. [31 An important contribution 
to the bound comes, of course, from the CHOOZ reactor experiment ^ combined with 
the determination of jAml^l from atmospheric and long-baseline experiments. Due to 
a complementarity of low and high energy solar neutrino data, as well as solar and 



Three-flavour neutrino oscillation update 



CD 

N 

c 
'(/) 



0.15 



P 0.1 



0.05 



1 1 1 
— f 


1 1 1 1 1 1 1 1 
—-^KamLAND 


1 1 1 1 


- 


\ global 


_ 




solar - 










M 




1 1 1 


1 wwrim, 


1 



Q I I I I I I <Vi \ I LL^m I I I 

0.1 0.2 0.3 0.4 0.5 



Sin 



12 



20 



15 



5 - 







~\ — I — I — r 



T — I — I — r 



- 2008 
-- 2007 



3a 



/ 



/ 

// 



// 



/ 
// 



~7n 



- 90% CL / 



0.05 



0.1 



Sin 



13 



Figure 4. Left: Allowed regions in the {612-613) plane at 90% and 99.73% CL (2 dof) 
for solar and KamLAND, as well as the 99.73% CL region for the combined analysis. 
Am|j^ is fixed at its best fit point. The dot, star, and diamond indicate the best fit 
points of solar, KamLAND, and combined data, respectively. Right: profile from 
solar and KamLAND data with and without the 2008 SNO-NCD data. 



KamLAND data, one finds that also solar+KamLAND provide a non-trivial constraint 
on ^13, see e.g., Refs. 0,3 Q- We obtain at 90% CL (3a) the following limits Q: 



0.060 (0.089) (solar+KamLAND) 
sin^^i3<<J 0.027 (0.058) (CH00Z+atm+K2K+MIN0S) (3) 
0.035 (0.056) (global data) 

In the global analysis we find a slight weakening of the upper bound on sin^ ^13 due to 



the new data from 0.04 (see Ref. 33| or v5 of [8|) to 0.056 at 3a. The reason for this is 
two-fold. First, the shift of the allowed range for |Am3]^| to lower values due to the new 
MINOS data implies a slightly weaker constraint on sin^ 613 (cf. Fig. |3l left), and second, 
the combination of solar and new KamLAND data prefers a slightly non-zero value of 
sin^ 6'i3 which, though not statistically significant, also results in a weaker constraint in 
the global fit (cf. Fig. El right). 



As has been noted in Ref. 16|] the slight downward shift of the SNO CC/NC ratio 
due to the SNO-NCD data leads to a "hint" for a non-zero value of ^13. From the 
combination of solar and KamLAND data we find a best fit value of sin^ 613 = 0.03 
with Ax^ = 2.2 for 6*13 = which corresponds to a 1.5a effect (86% CL). We illustrate 
the interplay of solar and KamLAND data in the left panel of Fig. HI The survival 
probability in KamLAND is given by 

pKamL ^ ^^^4 _ ^^^2 3^^^ gJ^S , (4) 

J Note that the bounds given in Eq. ([3|) are obtained for 1 dof, whereas in Fig. [3] (left) the 90% CL 
regions for 2 dof are shown. 
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leading to an anti-correlation of sin^ and sin^ 612 
for solar neutrinos one has 



see also IJ, [34 . In contrast 



psolar ^ 
ee ^ 



cos^On (1 - isin2 2^i2) 
cos^ 6^13 sin^ 6^12 



low energies 
high energies 



(5) 



Eq. ([5]) shows a similar anti-correlation as in KamLAND for the vacuum oscillations of 
low energy solar neutrinos. For the high energy part of the spectrum, which undergoes 
the adiabatic MSW conversion inside the sun and which is subject to the SNO CC/NC 
measurement, a positive correlation of sin^^is and sin^ ^12 emerges. As visible from 
Fig. m (left) and as discussed already in Refs. jsl, Is^], this complementarity leads to a 
non-trivial constraint on 613 and it allows to understand the hint for a non-zero value 
of ^13, which helps to reconcile the slightly different best fit points for 612 for solar and 
KamLAND separately , |l6|. This trend was visible already in pre-SNO-NCD data, 
though with a significance of only 0.8a, see Fig. H] (right) showing the present result 
together with our previous one from v6 of jsl. 

Let us briefly comment on a possible additional hint for a non-zero ^13 from 



atmospheric neutrino data 15|, l30|; Refs. 16|, [SOj find from atmospheric+long- 
baseline+CHOOZ data a 0.9a hint for a non-zero value: sin^6'i3 = 0.012 ±0.013. In our 
atmospheric neutrino analysis (neglecting Amgi effects) combined with CHOOZ data 



the best fit occurs for ^13 = (cf . Fig. [3l right), in agreement with Ref . [29| . Also, in 



the atmospheric neutrino analysis from Ref. SlJ (which does include Am|]^ effects, as 



Refs. 16|, ISOj) the preference for a non-zero 6*13 is much weaker than the one from [30 
with a Ax^ ^ 0.2. In our global analysis the hint from solar+KamLAND gets diluted 
by the constraint coming from atmospheric-|-CHOOZ data, and we find the global 
minimum at sin^^i3 = 0.01, but with ^13 = allowed at 0.9cr (A^^ = 0.87). Hence, we 
conclude that at present there is no significant hint for a non-zero ^13. As already stated, 
the origin of slightly different conclusions of other studies is related with including or 
neglecting the effect of solar terms in the atmospheric neutrino oscillation analysis, and 
translates also into a possibly nonmaximal best fit value for 623. Note, however, that 
all analyses agree within Ax^ values of order 1 and therefore there is no significant 
disagreement. A critical discussion of the impact of sub-leading effects in atmospheric 
data on Oi^ and 6*23 as well as a comparison of the results of different groups can be 



found in Ref. 33 



Before summarizing let us update also the determination of the ratio of the two 
mass- squared differences. 



a = 



_Am|i_ 

\Aml^\ 



0.032, 0.027 < a < 0.038 (3a) 



(6) 



which is relevant for the description of CP violation in neutrino oscillations in long- 
baseline experiments. 
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4. Summary 

In this work we have provided an update on the status of three-flavour neutrino 
oscillations, taking into account the latest available world neutrino oscillation data 
presented at the Neutrino 2008 Conference. Our results are summarized in Figures [11 
Eland El Table m provides best fit points, la errors, and the allowed intervals at 2 and 
3cr for the three-flavour oscillation parameters. 



parameter 


best fit 


2(7 


3(T 


Am2^ [lO-^eV^] 


' ■'-'^-0.20 


7.25-8.11 


7.05-8.34 


|Am|i| [lO-^^eV^] 


2 4n+°-^2 


2.18-2.64 


2.07-2.75 


sin^ 6^12 


n Qn4+°-°22 

U.OU^_0.oi6 


0.27-0.35 


0.25-0.37 


sin^ 6^23 




0.39-0.63 


0.36-0.67 


sin^ 6113 


U-Ui_o.011 


< 0.040 


< 0.056 



Table 1. Best-fit values with Icr errors, and 2a and 3ct intervals (1 d.o.f.) for 
the three-flavour neutrino oscillation parameters from global data including solar, 
atmospheric, reactor (KamLAND and CHOOZ) and accelerator (K2K and MINOS) 
experiments. 
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Appendix A. Updated analysis as of February 2010 

Appendix A.l. Updates in the solar neutrino analysis 



SSM : We consider the recently updated standard solar model from 35(]. Among the 



different models presented in that reference, we use the low metallicity model, labelled 



36 



as 



as AGSS09, that incorporates the most recent determination of solar abundances 
our standard choice. The solar abundances in that model are a bit higher than previous 
determinations by the same group, alleviating the disagreement with helioseismic data. 
From the point of view of solar neutrinos, the most important changes with respect to 



the previous SSM used in our analysis (BS05(OP), with high metallicities [37|) is the 
15% and 5% reduction in the Boron and Beryllium fluxes respectively. This is due to 
the reduced central temperature in the new model with respect to the previous one. 
Given the condition of fixed solar luminosity, this reduction is compensated by a slight 
increase in the pp and pep neutrino fluxes. We discuss also the impact of a new SSM 
with high metallicity, the GS98 model (presented in [ssf as well), see also the recent 



discussion in 38 . 
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SAGE : In our present analysis we have updated the capture rate of solar 



n. 



neutrinos measured in SAGE: 65.4;^3;g;*;2 s SNU [39], compared to the previous value 



ee.g^gllg s SNU |40|. Note that the recently published reanalysis of GALLEX data [41 



has been reported already at Neutrino 2008 13|] and was therefore included in the 
original version of this paper. 

SNO : In our update we include also the results from the recent joint re-analysis of data 
from the Phase I and Phase II (the pure D2O and salt phases) of the Sudbury Neutrino 
Observatory (SNO) 42]. In this analysis, an effective electron kinetic energy threshold 
of 3.5 MeV has been used (Low Energy Threshold Analysis, LETA), and the total flux 
of '^B neutrinos has been determined to be 

<j)Mc = 5.140l°:i^°(stat)+S:;??(syst) x 10^ cm'^s'^ . (A.l) 

Comparing this number with the result obtained in SNO j'hase III (the NC detector 
phase, NCD): 5.54l°J?(stat):^[];34(syst) x 10^ cm-^g-i [Iflj, one can see that the 
determination of the total neutrino flux has been improved by about a factor 2. These 
improvements have been possible thanks mainly to the increased statistics, in particular 
the NC event sample in the LETA is increased by about 70%, since the previously used 
higher energy thresholds of 5 MeV in phase I and 5.5 MeV in phase II have cut away 
a significant portion of the NC events. Furthermore, energy resolution, backgrounds 
suppression, and systematic uncertainties have been improved. We include the LETA 
SNO data by fitting the predicted energy- dependent neutrino survival probability and 
day-night asymmetry in terms of the polynomials given by the SNO collaboration, see 



Tabs. XXVI and XXVII in j42j. We have checked that our results agree with the analysis 
including all solar neutrino experiments made by SNO. Note that we have adopted in 
our present analysis of the SNO-NCD phase data the detailed correlations between the 



CC, NC and ES neutrino fluxes recently given by the SNO collaboration 431], thereby 
improving our previous treatment presented in Sec. 2. 

Fig. lAll shows the impact of the updates in the solar analysis on the determination 
of the solar parameters. The left panel compares solar and solar+KamLAND allowed 
regions for the previous and updated analyses, the middle panel illustrates the impact 
of the SNO LETA analysis, and the right panel shows the effect of changing between 
the low (AGSS09) and high (GS98) metallicity solar models. We observe that the main 
changes come from the SNO LETA analysis, whereas the impact of solar metallicity is 
small. In general changes are rather small, once KamLAND data is added to the solar 
data, with a small tightening of the lower bound on sin^^^i2. The best fit point value 
for the solar mixing angle has been shifted to a slightly higher value mainly due to 
the lower values reported for the total boron neutrino flux, either from the SNO LETA 
measurements as well as from the updated SSM with low metallicities. The allowed 
range for Am^i of solar only data has been somewhat reduced, however this effect gets 
completely diluted after combining with KamLAND data. The updated best fit values 
and allowed ranges for sin^ 612 and Am^i can be found in Tab. IA1[ 

The impact of the updated solar analysis on 613 is illustrated in Fig. IA2[ Our 
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Figure Al. Impact of the changes in the solar neutrino analysis. In all panels the 
blue (shaded) regions corresponds to the 3cr regions from solar and solar+KamLAND 
updated analysis. The regions delimited by the red contour curves correspond to our 
previous analysis (left), an analysis using the previous high-threshold SNO phase I 
and II analysis but the same solar model (middle), and an analysis using the high 
metallicity GS98 instead of our standard low metallicity AGSS09 solar model (right). 



analysis of solar + KamLAND data gives sin^^is = 0.022^Q Qjf in excellent agreement 



0.015 

with the value obtained by the SNO Collaboration [42[, sin^^ig = 0.0200+o:oi63- Hence, 
we obtain a lower best fit value with respect to the one we obtained in our previous 
analysis (sin^ ^13 = 0.03). This is due to the fact that now solar data prefer a somewhat 
higher value of 6*12 (as KamLAND does), and therefore, a smaller value of 613 is required 
to reconcile solar and KamLAND data, as can be seen by comparing left panels of 
Figs, m and IA2I The fact that now solar data prefer a larger value for sin^ 6*12 results 
in a stronger bound on ^13 from the combination of solar + KamLAND data. The 
allowed solar region in the panel (sin^ Ou, sin^ ^13) is more shifted to the right (because 
of the higher 612 preferred by the new smaller boron neutrino flux), where the allowed 
KamLAND region is narrower. At ^13 = we find Ax^ = 2.2, same value as before. 

As stated above, the small improvement in the ^13 bound is related to the solar 
model used. For models with higher solar metallicities like GS98, a slightly weaker 



bound is obtained 38|], see Fig. IA2I (right). In that case we obtain a slightly larger best 
fit point, sin2^i3 = 0.027l|j:[|}^ and Ax^ = 3.05 at ^13 = 0. 

Appendix A. 2. MINOS v^, appearance data 

In Ref. [4j] a search for transitions by the MINOS experiment has been 



presented, based on a 3.14 x 10^*^ protons-on-target exposure in the Fermilab NuMI 
beam. 35 events have been observed in the far detector with a background of 
27 ± 5(stat) ± 2(syst) events predicted by the measurements in the near detector. This 
corresponds to an excess of about 1.5cr which can be interpreted as a weak hint for 
Ue appearance due to a non-zero ^13. We fit the MINOS z/g spectrum by using the 



GLoBES simulation software |45|, where we calibrate our predicted spectrum by using 
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Figure A2. Left (update of Fig. |4j left): Allowed regions in the {9i2 — O13) plane 
at 90% and 99.73% CL (2 dot) for solar and KamLAND, as well as the 99.73% CL 
region for the combined analysis. Right: Ax^ as a function of sin^ 6*13 . The blue curves 
illustrate the impact of the updates in the solar neutrino analysis on the bound from 
the global solar+KamLAND data. The red curves show the constraint coming from 
the MINOS Ve appearance (red solid) and NC (red dashed) data, where we show the 
Ax^ assuming NH (thin) and IH (thick) , both with respect to the common minimum, 
which occurs for IH. The green solid curve corresponds the bound from CHOOZ + 
atmospheric + K2K + MINOS (disappearance) data. 



the information given in 46| . A full three-flavour fit is performed taking into account a 
7.3% uncertainty on the background normalization (Tab. I of i^l), and a 5% uncertainty 
on the matter density along the neutrino path. 

In the MINOS detector, being optimized for muons, it is rather difficult to identify 
z/g CC events, since they lead to an electromagnetic shower. NC and z/^ CC events often 
have a similar signature, and hence lead to a background for the Uf, appearance search. 
Indeed, in Ref. 47|] an analysis of "NC events" has been performed, where "NC events" 
in fact include also Uf, CC events due to the similar event topology. Therefore, a possible 
i^fj. i^e oscillation signal would contribute to the "NC event" sample of 47j and these 
data can be used to constrain 9i^. We have performed a fit to the observed spectrum, 
again using the GLoBES software, by summing the NC events induced from the total 
neutrino fiux with the Ue CC appearance signal due to oscillations. We include a 4% 
error on the predicted NC spectrum and a 3% error on the CC induced background 
(Tab. II of 0). 

In Fig. IA2I (right) we show the constraint on sin^ ^^13 from these MINOS data. The 
has been marginalized with respect to all parameters except ^13, where for the solar 
and atmospheric parameters we imposed Gaussian errors taken from Tab. without 
including any other information on ^13 except from MINOS. We show the Ax^ profiles 
for Ve appearance data and NC data, for a fixed neutrino mass hierarchy. The best fit 
point is always obtained for the inverted hierarchy (IH, ^m\^ < 0), and in that case 
in general the constraint on sin^ ^13 is weaker, since for IH the matter effect tends to 
suppress the v^, appearance probability. The Ax^ for normal hierarchy (NH, Amg^^ > 0) 
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Figure A3. The constraint on sin^ ^13 from MINOS appearance data, solar + 
KamLAND data, atmospheric + CHOOZ + K2K + MINOS (disappearance as well as 
appearance), and the combined global data. 

is given with respect to the best fit for IH. In the global analysis we also marginalize 
over the two hierarchies, and hence, the actual information from MINOS comes from 
the IH. 

We see from the figure that MINOS Uf, appearance data shows a slight preference 
for a non-zero value of ^13, with a best fit point of sin^ ^13 = 0.032(0.043) for NH (IH) 
with Ax^ = 1.8 at sin^ ^13 = 0. In contrast, no indication for a non-zero ^13 comes from 
the NC data. Furthermore, one observes that NC gives a slightly more constraining 
upper bound on sin^ ^13 than z/g appearance, while both are significantly weaker than 
the bound from i/^ disappearance data + CHOOZ or solar+KamLAND. Let us mention 
that the result for the NC analysis strongly depends on the value assumed for the 
systematic uncertainty, whereas the u^. appearance result is more robust with respect to 
systematics, being dominated by statistics. 

In the global analysis we do not combine the x^'s from MINOS z/g and NC data, 
since presumably the data are not independent and adding them would imply a double 
counting of the same data. Therefore, we adopt the conservative approach and use only 
z/g appearance data without the information from NC data in the global analysis. We 
have checked, however, that adding both MINOS data sets leads to practically the same 
result in the global fit, both for the "hint" for ^^13 > as well as the global bound, the 
latter being dominated by other data sets. 
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Appendix A. 3. Updated global three-flavour fit 



The present situation on the mixing angle ^13 is summarized in Fig. I A3 1 We obtain the 
following bounds at 90% {3a) CL: 



sin^ ^13 < 



0.053 (0.078) 
0.033 (0.061) 
0.034 (0.053) 



(solar+KamLAND) 

(CH00Z+atm+K2K+MIN0S) (A.2) 
(global data) 



We note a slight tightening of the bounds from solar+KamLAND as well as the global 



bound, due to the update in the solar analysis, see [Appendix A.l[ whereas the bound 
from CH00Z+atm+K2K+MIN0S gets slightly weaker, due to MINOS appearance 
data. In the global analysis we obtain the following best fit value and la range: 

sin2^i3 = 0.013^°:°^^ (A.3) 

This corresponds to a 1.5a hint for 6*13 > (Ax^ = 2.3 at 6^13 = 0). As discussed in sec. [3] 
above, in our previous analysis the 1.5cr hint for 613 > from solar+KamLAND data 
was diluted after the combination with atmospheric, long-baseline and CHOOZ data, 
resulting in a combined effect of 0.9a. Now, thanks to the new MINOS appearance 
data, we find that the atmospheric + long-baseline + CHOOZ analysis already gives a 
nonzero best fit value of 6*13 (see Fig. IA3I) . leading to the above global result, eq. IA.3[ 

Finally, let us comment on the possible hint for a non-zero ^13 from atmospheric 
data [la, l30|, as discussed in sec. [3l The possible origin of such a hint has been 



investigated in Ref. 



48 



and recently in 



38 



see also 49 . From these results one 



may conclude that the statistical relevance of the claimed hint for non-zero ^13 from 
atmospheric data depends strongly on the details of the rate calculations and of the 
analysis. Furthermore, the origin of that effect might be traced back to a small excess 
(at the la level) in the multi-GeV e-like data sample in SK-I, which however, is no 
longer present in the combined SK-I and SK-II, as well as SK-I+II+III data. 



parameter 


best fit 


2a 


3a 


Aml^ [lO-^eV^] 




7.22-8.03 


7.03-8.27 


lAml^l [lO-^^eV^] 


2 40+°-^2 


2.18-2.64 


2.07-2.75 


sin^ 612 


n c+0.019 


0.29-0.36 


0.27-0.38 


sin^ 6^23 




0.39-0.63 


0.36-0.67 


sin^ 6^13 


Q Q-|O+0.013 

U.UiO_o.o09 


< 0.039 


< 0.053 



Table Al. Current update of Tab. [TJ Best-fit values with la errors, and 2a and 3a 
intervals (f d.o.f.) for the three-flavour neutrino oscillation parameters from global 
data including solar, atmospheric, reactor (KamLAND and CHOOZ) and accelerator 
(K2K and MINOS) experiments. 



Tab. lAll gives an updated summary of the present best fit values and allowed ranges 
for the three-flavor oscillation parameters. 
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